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FAST MULTTPOLE METHOD SOLUTION 
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University of Illinois 
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1. Introduction 

The fast multipole method (FMM) [1-6] speeds up the matrix- vector mul- 
tiply in the conjugate gradient (CG) method when it is used to solve the 
matrix equation iteratively. In this paper, FMM is applied to solve the elec- 
tromagnetic scattering from three dimensional arbitrary shape conducting 
bodies. The electric field integral equation (EFIE), magnetic field integral 
equation (MFLE), and combined field integral equation (CFIE) are consid- 
ered. FMM formula for CFIE has been derived, which reduces the complexity 
of a matrix-vector multiply from 0(N 2 ) to ©(AT 1 - 5 ), where N is the number 
of unknowns. With a nonnested method, using the ray-propagation fast mul- 
tipole algorithm (RPFMA), the cost of a FMM matrix-vector multiply is 
reduced to 0(JV 4 ' 3 ). We have implemented a multilevel fast multipole algo- 
rithm (MLFMA), whose complexity is further reduced to O(MogJV). The 
FMM also requires less memory, and hence, can solve a larger problem on a 
small computer. 

2. The Fast Multipole Method (FMM) 

Practical electromagnetic problems are often three-dimensional and in- 
volve arbitrary geometry. The arbitrary surface is described by dividing it 
into a number of connected patches which are mathematically described as 
parametric quadratic surfaces [7]. For conducting objects, the electric field 
integral equation (EFIE) is given by 

< • J G(r, r') • JOW = • E'"(r), (1) 

and magnetic field integral equation (MFIE) for closed conducting objects is 
given by 

27rt - J(r) - i • rt x V x J dS'g(ry)3(r') = 4jrt -ft x IT(r), (2) 

where 

G(r, rO = (I - j^W'Mr, «0, g(r,r') = (3) 

t This work was supported by NASA under grant NASA NAG 2-871, Office of Naval Re- 
search under grant N00014-89-J1286, the Army Research Office under contract DAAL03- 
91-G-0339, and the National Science fbundation under grant NSf ECS 92-244C6. 
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Then combined field integral equation (CFIE) for closed conducting objects 
is simply a linear combination of EFIE and MFIE. The integral equations 
are approximated by matrix equations using the method of moments (MOM) 
with specially designed basis functions for subdomains which contain surface < 
curvature. The basis functions are the generalized rooftop functions. \ 

T he FMM idea is first to divide the subscatterers into groups. Then, \ 
addition theorem is used to translate the scattered field of different scattering \ 
centgjH-git hin a group iutu a single cental *?- Hence, the number of scattering / 
'centers is reduced. Si milarly, for each group, the field scattered by all th e / 
other group c enters can be first "receiwd" by tEe gro5p~center, -»""> *h°r / 
"redistributed" to the suhscatterers belonging to the group. / 
Aft er nemo d e rivati ons [2,6] using the addition theorem, we can rewrite 
the matrix- vector multiplication as 





m , (B m 



•eo ra 



Omm-(r mm . • k) = >T «'(2' + l)h?\kr mm .)P,(r mm . ■ k) 



V /m ,(fc) = a dSe*T*»(I - kk) ■ t f (r^) 

-(1 -a)k x j s dS^""t j {T jm ) x h 



" i(i) = i 




«iS e *""-'(I-fcJt)-j,(r im .) 



(7) 



he first te rm in f4) is the contribution from nearby groups (including the/ 

' "and the second-te rm is the far interaction calculated by r'MM. 
Secomputation cost using (4) with 2-level FMM is in order 0(1 
T o implemen t multilggel fust imiltipnlp plg™-i>k™ ( MTrFMAI. the entire 
object is first enclosed into a large cube, which is partitioned into eight smaller 
cubea. Ea ch subcube is then recursi vely subdivided' into 8nKdfe r~rii fres until - 

tKe edge length of the finest, ruhp is about half^if a w="" ,,tJ "S t \ When the 
cube Becomes larger trom the finest levej to thecoarsest level, the numbers of 
multipole expansions should increasefln t he first sweep, the outer multipole 
expansions are computed at the nneBTTlevel, Uieu tne expansions for larger 
cube are obtained using shifting and interpolation. At the coarsest level, the 
local multipole expansions contributed from well-separated cubes are calcu- 
lated using the second part of (4). At the second sweep, the local expansions 
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for smaller cube include the contributions from parent cube using, shifting 
and anterpolation, and from the well-separated cube at this level but not 
well-separated at the parent level. At the finest level, the contributions from 
non-wcll-separated cube are calculated directly. Since only nonempty cubes 
are considered, the complexity for MLFMA is further reduced to O(JVlogJV). 

3. Results and Conclusions 

Figure 1 shows the validation of the numerical result from combined field 
integral equation (CFIE) with FMM against the Mie series solution of the 
bistatic RCS of a metallic sphere of radius lm at frequency of 0.72GHz for 
the parallel polarization. 9408 unknowns with 2-levcl FMM are used. The 
solutions of CFEE with FMM agree with Mie series very well. 

Figure 2 shows the bistatic RCS of a one meter long metallic square plate 
at 4.5GHz in the xy plane with incident angle 8 = 45". 32512 unknowns with 
6-level FMM are used. The calculation is done by solving EFIE on a SUN- 
SPARC-2 with 64MB RAM. There is a good agreement between oiir results 
and the approximation by physical optics when the RCS is bigger than 0 dB. 

In conclusion, the fast multipolc method (FMM) has been implemented 
to speed up the matrix-vector multiply in the CG method when it is used to 
solve EFIE, MFIE, and CFIE. At all frequencies, CFIE has an unique solu- 
tion, and converges faster than EFIE and MFIE since the matrix from CFIE 
has a smaller condition number than those from EFIE and MFIE. FMM 
approach reduces the complexity of a matrix- vector multiply from 0(N 2 ) to 
0(N 1S ). With a multilevel fast multipole algorithm (MLFMA), the complex- 
ity is further reduced to O(NlogN). The FMM also requires less memory, 
and hence, can solve a larger problem on a small computer. 
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Figure 1. Validations of CFIE with FMM against the Mie 
series of the bistatic RCS of a metallic sphere of radius lm at 
0.72GHz for W polarization. The RCS is normalized by ira 2 . 

W=15X square plate, N=32512. 45° incidence, 6 levels 
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Figure 2. Bistatic RCS of a metallic square plate of length 
lm at 4.5GHz for W polarization with 45° incident angle. 
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